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We investigate the relative size of various twist-3 quark-gluon correlation functions relevant to 
single transverse spin asymmetries (SSAs) in a quark-diquark model of the nucleon. We calculate 
the quark-gluon correlation function T qi p(x, x) that is responsible for the gluonic pole contribution 
to the SSAs, as well as T q ,F(0,x) and Ta 9 ,f(0, x) responsible for the fermionic pole contributions. 
We find in both cases of a scalar diquark and an axial-vector diquark that at the first nontrivial 
order only the T qi F(x, x) is finite while all other quark-gluon correlation functions vanish. Using the 
same model, we evaluate quark Sivers function and discuss its relation to the T qi F{x,x). We also 
04 ■ discuss the implication of our finding to the phenomenological studies of the SSAs. 
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I. INTRODUCTION 

O The novel phenomenon of single transverse-spin asymmetry (SSA), An = (<r(S±) — (r(—S±))/((r(S±) + a(—S±)), 

defined as the ratio of the difference and the sum of the cross sections when the single transverse spin vector Sj_ is 
flipped, was first observed in the hadronic A production at Fermilab in 1976 as a surprise [l[. Large SSAs, as large as 
30 percent, have been consistently observed in various experiments involving one polarized hadron at different collision 
energies, and have attracted tremendous interest from both experimental and theoretical sides in recent years 0. The 
I \ size of the observed SSAs presented a challenge to the early QCD calculation *3j. As a consequence of the parity 
i and time-reversal invariance of the strong interaction, the SSAs in high energy collisions are directly connected to the 
transverse motion of quarks and gluons inside the transversely polarized hadron. The measurement of SSAs provides 
an excellent opportunity to probe a new domain of QCD dynamics. The understanding of the physics behind the 
measured asymmetries should have a profound impact on our knowledge of strong interaction and hadron structure. 
Two complementary QCD-based approaches have been proposed to analyze and to explore the physics behind 
t^J" ■ the measured SSAs: the transverse momentum dependent (TMD) factorization approach [j-tloj and the collinear 
factorization approach [lll - [l5j . Both approaches have been applied extensively to phenomenological studies (lfj| - [23 |. 
The TMD factorization approach is more suitable for evaluating the SSAs of scattering processes with two very 
different momentum transfers, Qi ^ Q2 ^ Aqcd- Having one large observed scale, Qi 3> Aqcd, is necessary for 
J> . using perturbative QCD and the TMD factorization approach although it is not sufficient [25l - [28j . For observables for 
which the TMD factorization is valid, this approach has an advantage for directly probing active parton's transverse 
motion at the scale, 0(Q-2), inside a polarized hadron in the form of TMD parton distribution functions (PDFs). 
On the other hand, the collinear factorization approach is more relevant to the SSAs of scattering cross sections 
with all observed momentum transfers Q 3> Aqcd- In the QCD collinear factorization approach, the leading power 
contribution to the cross sections in the 1 /Q expansion cancels in evaluating the asymmetry because of the parity and 
time-reversal invariance of the theory. Therefore, the asymmetry directly probes the correlation of quarks and gluons 
inside a polarized nucleon in the form of the twist-3 quark-gluon and tri-gluon correlation functions [l2| . Although the 
two approaches each have their own kinematic domain of validity, they describe the same physics and are consistent 
with each other in the regime where they both apply [3(| ■ 

In both TMD and collinear factorization approaches of QCD, the size of calculated SSAs is proportional to some 
nonperturbative functions: the TMD PDFs and the twist-3 three-parton correlation functions, respectively. The 
predictive power of both approaches relies on the validity of the respective factorization and the knowledge of these 
nonperturbative functions |3Xl , |32| . QCD perturbation theory could be used to study the quantum evolution of 



, — 1 

X 

5-H 



'Electronic address: zkang@bnl.gov 

^Electronic address: jqiu@bnl.gov 

t Electronic address: hongzh87@iastate.edu 



2 



these functions from one perturbative scale to another where these functions were probed ^33-36]. But, the absolute 
normalizations of these functions or the boundary conditions - the input functions for solving the evolution equations 
have to be extracted from data of measured asymmetries. With the recent measurements of SSAs [3714391 ] . we have 
gained valuable information on the TMD PDFs |40j and twist-3 correlation functions [22[ . Although precise data from 
future experiments could certainly help fix these nonperturbative functions, a good model calculation of these unknown 
functions could provide important insight into the mechanisms for generating the observed novel asymmetries and 
valuable guideline to the relative importance and size of various functions. In this paper, we present our calculations 
of all twist-3 quark-gluon correlation functions relevant to the SSAs in the collinear factorization approach in a quark- 
diquark model of the nucleon [f| |4l|. We calculate these quark-gluon correlation functions with both scalar and 
axial- vector diquarks. We also evaluate in the same model the quark Sivers function, the spin dependent part of the 
TMD quark distribution, and discuss its relation to the twist-3 quark-gluon correlation function, T q p(x,x) [9j. 

In order to generate a nonvanishing SSA in high energy hadronic collisions, one needs to generate a parton-level 
spin flip and a phase difference between the scattering amplitude and its complex conjugate. In the QCD collinear 
factorization approach to the SSAs, the spin flip at the hard collision was achieved by the interference between an active 
single parton state and an active two-parton composite state of the scattering amplitude; and the phase difference was 
generated by the interference between the real part and the imaginary part of the short-distance partonic scattering 
amplitude 12]. We obtain the leading contribution to the imaginary part of the partonic scattering amplitude by 
taking the unpinched pole of the partonic scattering amplitude [111 . Il2j . It is the interference between the single active 
parton state and the two-parton composite state that requires the calculated SSAs to be proportional to the twist-3 
quark-gluon correction functions, T qi p(x\, x 2 ) and T^ Qi p(xi 1 X2), and tri-gluon correlation functions, Tq'p\xi, x 2 ) 

and T^q F (xx, x 2 ), convoluted with corresponding partonic scattering through two independent momentum fractions 
of the three active partons, x\ and x 2 [HI, |36|. Taking the pole of the parton scattering amplitude effectively fixes 
one of the two momentum fractions. Depending on the number of observed hard momentum scales, the partonic 
scattering amplitude has different pole structure. For cross sections with a single observed hard scale, such as pr 
of single inclusive pion production in hadronic collisions, the leading pole contribution is from taking the residue of 
the pole, which is effectively setting the momentum fraction of one of the three active partons to zero [ill, El ■ This 
contribution is often referred as the soft-pole contribution. The so-called gluonic (or fermionic) pole contribution 
refers to the situation when the active gluon (or (anti)quark) momentum fraction was set to zero. For cross sections 
with more than one observed hard scale, such as inclusive pion production in lepton-hadron deep inelastic scattering 
when both pion momentum pt and virtual photon invariant mass Q are large, the leading pole contribution could also 
come from the situation when all active parton momentum fractions are finite, known as the hard-pole contribution 
[42I [43| . In this paper, we present the model calculation only for quark-gluon correlation functions corresponding to 
the soft-pole contribution to the SSAs, which include T q ^p(x, x) and T& qt p(x, x) for the gluonic pole contribution, and 
T qt p(0,x), T^. qt p(0,x), T Qt p(x,0), and Ta 9! _f(cc,0) for the fermionic pole contribution. 

In general, the calculated SSAs of cross sections with one observed hard momentum scale receive contributions 
from both t he g luonic and fermionic poles, and the partonic hard parts for these two contributions often have the 



similar size [22], |24|. On the other hand, the quark-gluon correlation functions corresponding to the gluonic pole 



and fermionic pole contribution represent very different dynamical structure inside the polarized proton. For the 
gluonic pole contribution, the quark-gluon correlation functions, T qi p(x,x) and T^ q ^p(x,x), represent the quantum 
interference between a quark state of momentum fraction x and a quark-gluon composite state of the same momentum 
fraction with the quark carrying all the momentum fraction x, while for the fermionic pole contribution, the total 
momentum fraction of the quark-gluon composite state is carried by the gluon. The relative size of these two types of 
quark-gluon correlation functions certainly provides interesting information on the dynamical structure of a polarized 
proton. 

In terms of the simple quark-diquark model of the nucleon 0, l4lj . we find in both cases of a scalar diquark 
and an axial-vector diquark that at the first non-trivial order, all quark-gluon correlation functions corresponding 
to the fermionic pole contribution, T Qt p (0,x), T/\ qt p(0, x), T qt p(x,0), and T& qi p(x, 0), vanish. For those functions 
corresponding to the gluonic pole contribution, T qt p(x,x) is finite while T^ q ^p{x 1 x) = which is consistent with the 
result of symmetry argument [33j]. Our results, although from a model calculation, indicate that the fermionic pole 
contribution to the SSAs is likely to be less important than the gluonic pole contribution. This conclusion seems to 
be consistent with a general expectation that a quark-gluon state with the quark carrying all of its momentum is 
more likely than a quark-gluon state with the gluon carrying all of its momentum to interfere with a quark state of 
the same momentum [12j . Our finding could help streamline the phenomenological studies of the SSAs by starting 
with a much smaller number of nonperturbative twist-3 correlation functions. 

The rest of our paper is organized as follows. In Sec. [Til we introduce the operator definition of all twist-3 quark- 
gluon correlation functions and discuss their symmetry properties. In Sec. Mil we introduce the quark-diquark model 
of the nucleon and its Feynman rules for both cases of a scalar diquark and an axial- vector diquark, and present our 
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calculations for the twist-3 quark-gluon correlation functions relevant to both gluonic and fermionic pole contributions 
to the SSAs. In terms of the same nucleon model, we calculate the quark Sivers functions in Sec. IIV1 and discuss 
the connection between the twist-3 quark-gluon correlation functions and the TMD parton distribution functions. 
Finally, we give our summary and conclusions in Sec. |V] 



II. THE TWIST-3 QUARK-GLUON CORRELATION FUNCTIONS 



The twist-3 three-parton correlation functions in the QCD collinear factorization approach to the SSAs could be 
represented by the cut forward scattering diagram in Fig. Q] with proper cut vertices [33J. These correlation functions 
measure the net effect of quantum interference between two scattering amplitudes of the transversely polarized proton: 
one with a single active parton and the other with two active partons, participating in the short-distance hard 
scattering. A complete set of twist-3 three-parton correlation functions relevant to the SSAs has been constructed 
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FIG. 1: The Feynman diagram representation for the twist-3 quark-gluon correlation functions, where ki « Xip with i = 1,2. 
Contracting the three active partons with different cut vertices leads to different quark-gluon correlation functions [331 ] - 

in Refs. [13, HH, which includes two independent quark-gluon correlation functions, T qi p(xi, x 2 ) and Ta<i.f(xi ,x 2 ). 
They could be derived from the following quark-gluon matrix element of a transversely polarized hadron (l4j , 

M°{x u x 2 ,s T ) = J ^^e-^ + ^ +l ^-^)P + ^(p, ST |^(0) 3 F CT +( 2/2 -)^(yr)|p, S T) 

= ~ [fie aSTnfl T q , F (xi,x 2 )+ 1 5 fii S ° T T Aq , F (x 1 ,x 2 ) + ---] , (1) 

where the proper gauge links have been suppressed and x% = ki ■ n/p ■ n with i = 1,2 are two independent parton 
momentum fractions, — [1 + ,0 - ,0jJ and = [0 + ,l~,0j_] are two light-like vectors with n-n = 1, and the 
ellipsis represents terms at twist-four and higher. From Eq. ([1]), we derive the expressions for the relevant quark-gluon 
correlation functions |33j |. 

T 9 , F (W) = / (2) 

w(zi,z 2 ) = I */^ e „ 1P+ , r+l( , 2 ^^ (3) 

From parity and time-reversal invariance, these two functions have the following symmetry property under the ex- 
change of the two arguments X\ <-> x 2 [H, HI]: 

Tg tF (xi,x 2 ) = Tg tF (x 2 ,xi), T A g >F (x 1 ,x 2 ) = -T AqiF (x 2 ,xi). (4) 

The leading order gluonic pole contribution to the SSAs is connected to the diagonal quark-gluon correlation func- 
tions, T qiF (x, x) and T Aq ^ F (x,x) (Til [l2T|. On the other hand, the leading order fermionic pole contribution to the 
SSAs is connected to the off-diagonal quark-gluon correlation functions, T qi p(0,x) and Ta 9i _f(0, x), or T qtF (x,0) and 
TAq.F(x,0) [TTl [T2|. From Eq. we have T Aq , F (x, x) = 0. In the next section, we calculate these correlation 
functions in the quark-diquark model of the nucleon, and test the symmetry properties in Eq. 



III. MODEL CALCULATION OF TWIST-3 QUARK-GLUON CORRELATION FUNCTIONS 

In this section, we calculate the twist-3 quark-gluon correlation functions relevant to the gluonic and fermionic pole 
contributions to the SSAs in the quark-diquark model of the nucleon |4l| . We consider two possible situations in 
which the spectator diquark is either a scalar particle or an axial-vector particle. 
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A. The quark-diquark model of the nucleon 

In the quark-diquark model of the nucleon |4l[ , the nucleon of mass M consists of a constituent quark of mass 
m and a diquark spectator of mass M s . The interaction between the nucleon, the quark, and the diquark is given by 





k 

(c) 

FIG. 2: Feynman diagrams to define the Feynman rules in the quark-diquark model of the nucleon: (a) vertex links the nucleon, 
the quark, and the diquark, (b) interaction vertex between the gluon and the diquark, and (c) the diquark propagator. The 
diquark could be a scalar particle or an axial-vector particle. The Lorentz indices are for the gluon and the axial-vector diquark. 

the following Feynman rule for the vertex in Fig. [2a), 

i\ s F s (k 2 ) scalar diquark, (5) 

iA|-y^ 5 F v (k 2 ) axial-vector diquark, (6) 

where A Sj „ represent the point-like interaction strength with subscripts s and v for a scalar and an axial-vector 
diquark, respectively, F s>v (k 2 ) are suitable form factors as a function of k 2 - invariant mass square of the constituent 
quark. F SiV (k 2 ) = 1 is for a point-like vertex interaction. As explained later, a properly chosen form factor could 
help control the ultraviolet behavior of the calculated quark-gluon correlation functions. The Feynman rule for the 
coupling between the gluon and the diquark in Fig. Hfb) is given by 

% g s (2p — 2k — q) T scalar diquark, (7) 

i g v V T ' ya (q,p — k — q,k — p) axial- vector diquark, (8) 

with the coupling strength g s and g v for a scalar and an axial- vector diquark, respectively. Here V Tla {q 1 p — k — q, k—p) 
is given by [4l| 

V Tja (q,p — k-q,k-p) = g T ~<(2q - p + k) a + g^ a (2p -2k- q) T + g aT (k - p - q)* 1 . (9) 

The Feynman rule for a scalar diquark propagator is the same as that of a normal scalar particle, while the Feynman 
rule for the axial- vector diquark propagator in Fig. [He) is given by 

d a0 (k,M s ), (10) 



k 2 - M 2 

where the polarization tensor d a @ (k, M s ) has the following form [4l[ 

n ■ k (n ■ ky 

which has the property, k a d a @(k, M s ) = when k 2 = M 2 . 

As we will show below, the twist-3 quark-gluon correlation functions calculated with the point-like coupling between 
the nucleon, the quark, and the spectator diquark, F SiV (k 2 ) = 1, have logarithmic ultraviolet divergences when 
k 2 — > oo. Since we are mainly interested in the long-distance behavior of the quark-gluon correlation functions, we 
could choose a proper form factor to eliminate the divergence from the region of phase space where k 2 3> M 2 , the mass 
scale of the nucleon, while preserve the dynamics at k 2 ~ M 2 . Several choices for the form factor were introduced 
and discussed in Ref. [4l|. In our calculation below, we assume that the form factor for a scalar diquark is the same 
as that for an axial-vector diquark, and choose a dipolar form factor [4l| 

F(k 2 ) = F s (k 2 ) = F v (k 2 ) = ^ " A 2 S , (12) 

where A^ > M 2 is an ultraviolet cutoff. Note that in Eq. (TT2l we multiplied the dipolar form factor in Ref. [4JJ by an 
extra so that the form factor has the right dimension. Such a difference by a constant factor does not affect any 
of our conclusions derived below. We will also demonstrate below that the introduction of this form factor smoothly 
suppresses the influence of the ultraviolet region of k\ or k 2 without affecting the main conclusions of this paper. 
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B. Calculation with a scalar diquark 



All quark-gluon correlation functions could be represented by the same cut forward scattering diagram in Fig. [T] 
The difference of various quark-gluon correlation functions is from the difference in cut vertices contracted to the 
three active partons in the diagram [33[. The form of cut vertices for both T q p(xi . x?) and T& q ^p(xi, x 2 ) as well as 
the tri-gluon correlation functions is gauge dependent and was derived in Ref. |33j . In this paper, we work in the 
light-cone gauge. For X\ = x + y and x 2 = x, the cut vertices are given by [33] 

v£5 = £ w (« - k -f) v s (y $) ^ nn ) [-9^] c q , (is) 

VaIf = ^ 2*gS (x-^jy6(y-^j h<wl C q (14) 

where g with g 2 = 47ra s is the strong coupling constant included in the definition in Eq. ([T]), C q is the fermionic color 
contraction factor given by [33j 

(C q f l3 = (t% , (15) 

with quark and gluon color indices, i,j = 1,2,3 = iV c and c = I, 2,..., 8 = N 2 — 1, respectively, and t c are the 
generators of the fundamental representation of SU (3) color. 




FIG. 3: The lowest order Feynman diagram for twist-3 quark-gluon correlation functions in the quark-diquark model. 



The contribution to the twist-3 quark-gluon correlation functions T q p{x + y, x) and T& q p(x + y, x) at the lowest 
non-trivial order is given by the Feynman diagram in Fig. [3] We first study these correlation functions with a scalar 
diquark. Applying the cut vertex in Eq. (JT5J) to the diagram in Fig. [3J we obtain 



+ m)] 

-F(k 2 )F((k + qf), (16) 



xe STanfl (2p -2k- q) T d aT {q) Tr [ 7 +(jt + i + m)(#+ M)j 5 fr(ft + m) 
III I 

V 

k 2 — m 2 — ie q 2 + ie (k + q) 2 — m 2 + ie (p — k — q) 2 

where e = + represents a small positive parameter, the superscript (s) indicates the scalar diquark, and the gluon 
polarization tensor d aT (q) is given by 

dar\q) = -gar H ■ (17) 

q ■ n 

Performing the integration over fc + , k~ , and q + by using the three (^-functions in Eq. (j!6p . we obtain 

T {s) ( , s _ _ M r gA * 9s f d2q± d2k± 1 f dq " 

i q , F {x + y,x) - "c^F mnp+ J (2^)2(2^)2^.2 +jL 2( m2)< / 2n 

xe ST°nn {2p _ 2k _ q y( q + gaT _ qa n T )Tr [ 7 +(jg + ^ + m)^ + M)j 5 Mft + m)] 

-F(k 2 )F((k + q) 2 ), (18) 



q 2 + ie (k + q) 2 — to 2 + ie (p — k — q) 2 — M 2 + i 

where 



in 



2 ^ ["1,2 , t2/2 



1-x 



[kl+L 2 (m 2 )} (19) 
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with L 2 s (m 2 ) given by 



L 2 s {m 2 ) = xM 2 + (1 - x)m 2 - x(l - x)M 2 



(20) 



independent of q~ . 

The integration over q~ is crucial and is done by taking the residue of relevant pole(s) of the integrand in Eq. (|18l) . 
which provides the necessary phase for a real quark-gluon correlation function T q .p(x + y, x). Since we are interested 
in the leading gluonic and fermionic pole contribution to the SSAs, we exam below the pole structure of the integrand 
in Eq. (|T5|) at y = (gluonic pole) and x + y = (fermionic pole) while x > 0. From 



(p-k- qf 



Mi 



-2(1 



y)p + q~ 



y(k 2 A 



M 2 



1-x 



2k± ■ q± 



ie 



0. 



and x + y < 1, we derive the location of the corresponding pole at 



1 



q = 



2(1 



y)p^ 



y(k 2 ± + M 2 ) 



2k± ■ qj 



(21) 



(22) 



which is in the upper half plane of the q . This pole survives and stays in the upper half plane at both limits: y = 



(gluonic pole) and x + y — (fermionic pole). However, the potential poles from q 



and (k + q) 2 -m 2 + ie = 



are sensitive to these two limits. For the quark-gluon correlation functions relevant to the leading fermionic pole 
contribution to the SSAs, we consider the pole structure at x + y = while y < since x > and find that 



2yp + q 



9a 







(23) 



provides a pole at 



2\y\pl 



in the upper half plane of the q , while 

(k + q) 2 - m 2 + ie = 2{x + y)p + (k + q)- - (k± + q±) 2 - m 2 + ie = -(k ± + q^) 2 - m 2 



(24) 



(25) 



does not contribute to any pole in the ^'-integration. That is, when x + y = and x > 0, the integrand of q~- 
intcgration in Eq. (fl8|) has two poles from (p — k — q) 2 — M 2 + ie = and q 2 + ie = and both of them are in 
the upper half plane of q~ . Since the integration dq~ in Eq. (|18|) is sufficiently converging when \q~\ — > oo, the q~ 

integration vanishes by closing the g _ -contour through the lower half plane. In conclusion, T^ S p{0,x) — from this 

leading order calculation with a scalar diquark, so as T^ F (x,0) = 0, which can be derived by an explicit calculation 
or the symmetry property T qt p(x,0) = T qt p(0,x). 

Now we turn to the limit at y = 0, which is relevant to the leading gluonic pole contribution to the SSAs. At y = 0, 
the pole structure of the q~ -integration in Eq. (TT8l) changes. At y = 0, 



(k + qf 
leads to a pole at 



m + ie = 2xp + q 



[(1 - x)M 2 -k\- M 2 } - (k 2 ± + qi _f 



m, 



0, 



q = 



2^+ 



1 



in the lower half plane of q , while 



2(1 - x)p+ 



ie = -q A 



[(1 - x)M 2 -k\- M 2 } - ie 



(26) 



(27) 



(28) 



is independent of q . Therefore, for the quark-gluon correlation functions relevant to the leading gluonic pole con- 



tribution to the SSAs, the integration of dq in Eq. (TT51) has two poles from (p — k — q)* 



M 



r2 



ie = and 



(fc + q) 2 —mr+i 

upper or the lower half plane, we obtain 
N c C F gX 2 s g 



with one in upper and one in lower half plane of q . By closing the q -contour in either the 

d 2 k 1 _ d 2 q ± [q\ - (q ± ■ S± ) 2 ] F{k 2 ) F((k + q) 2 ) 



4tt 



(1 - x)(m + xM) 



(2tt) 2 (2tt) 2 q\ [k\ + L 2 (m 2 )} [{k± + q±) 2 + L 2 (m 2 )} ' 



(29) 
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where k 2 is given in Eq. (|19[) and 



(k + q) 2 = 



[( k±+qx ) 2 +L 2 (m 2 )] 



(30) 



with L 2 s (m 2 ) given in Eq. (l20l) . 

The integration over the transverse momenta in Eq. (|29[) depends on the choice of the form factor. If we set 
F(k 2 ) = F((k + q) 2 ) = 1 for the point-like interaction between the nucleon, the constituent quark, and the spectator 
diquark, we obtain after integrating over d 2 k±, 



point— like 



N c C F g\ 2 g s 
16tt 2 

N c C F g\ 2 s9s 
16tt 2 



(1 - x){m + xM) 



(1 - x)(m + xM) 



da 
d 2 qi_ 



d 2 q± 



Q± - (U ■ s ±-) 2 



(2tt) 2 q 2 ± [a(l 



a)ql+L 2 (m 2 )} 



lnV^^S±^, (3 l) 



(2tt) 2 q± ^ql+AL 2 (m 2 ) ^qi+AL 2 (m 2 ) 



q± 



which has the logarithmic ultraviolet divergence from the region \q± \ — > oo. Since we are interested in the dynamics 
at the hadronic scale, we could use the dipolar form factor in Eq. (|12[) to remove the ultraviolet divergence. Using 
Eqs. (QU and ([3"0]). we have 



F(k 2 )F((k + q) 2 ) = (l-x) 2 (A 2 s ) 



2 fc 2 +L 2 (m 2 ) (k ± +q ± ) 2 + L 2 (m 2 ) 
[kl+L 2 (A 2 )] 2 [(k ± + q ± ) 2 + L 2 (A 2 )) 27 



(32) 



thus from Eq. ([29 



dipolar 0(Z7Tj^ 

= N c C F g\ 2 g s 
16(2tt) 3 



(1 - cc) 3 (m + xM) 



( A 2 



(2vr) 2 [fc 2 +L 2 (A 2 )] 



V L KA 2 



(33) 



where L 2 {A 2 ) is given in Eq. (|20p with m 2 replaced by the cutoff scale A 2 . Note that the form factor in Eq. (|3"2")l 
suppresses the ultraviolet contribution to the integration in Eq. (j29|) without altering the pole structure of the original 



diagram. Therefore, our general conclusion on Tff F (0,x) — T^ S p{x,Q) = remains whether we use the dipolar form 
factor or not. 

The calculation for T^ q ^p(x + y, x) is identical to that of T q ^p(x + y, x) except the cut vertex is replaced by the one 
in Eq. (| 14|) . Since the pole structure of the diagram is exactly the same, we obtain the same result for correlation 
functions relevant to the fermionic pole contribution, 



T^l F (Q,x) = -T^ F (x,0) = 



0) 



(34) 



From the symmetry property in Eq. we have expected the diagonal correlation function F (x,x) relevant to 
the gluonic pole contribution to vanish. As a consistent test of our model calculation, we verify this result explicitly 
as follows. Following the same procedure used to evaluate T^ s F {x,x) above, we use first the 5- functions to integrate 
over k + , k~ , q + , then the pole structure to integrate over q~ to get the necessary phase, and we obtain 



T 



Aq,F 



(x, x) 



N c C F g\ 2 s g s 

47T 



(1 -x)(m + xM) 



d 2 k 1 _ d 2 q ± q ± ■ s ± [2fc ± ■ s± + g ± ■ Sl _]F{k 2 ) F((k + qf 
(2^) 2 (2tt) 2 q\ \k\ + L 2 (m 2 )][(k ± + q ± ) 2 + L 2 {m 2 )] 



(35) 



We first consider the point-like interaction case setting F(k 2 )F((k + q) 2 ) = 1. Using the Feynman parametrization 
to combine the k±_ dependent denominator, we obtain 



Aq,F 



(x,x) 



N c C F g\ 2 s g s 



point— like 



4tt 



(1 - x)(m + xM) 



d 2 q ± 



da- 



(l-2a)(q ± -s ± ) 2 



(2tt) 2 (2tt) 2 J q l [jp ± + a (l - a)q 2 ± + L 2 (m 2 )Y 



= 0. 



(36) 



Here, the second line is due to the fact that the numerator of the a integral is antisymmetric under a 1 — a while 
the denominator and the integration limits are symmetric. From Eq. (|3"2"j) . it is clear that the inclusion of the dipolar 
form factor does not change the main feature of the a-dependence of the combined denominator, 



T 



Aq,F 



(x,x) 



dipolar 



N c C F g\ 2 s g s 3 2 x2 
— (1 - xY(m + xM)(A s ) 



d 2 q ± d 2 



da 



3!a(l-a)(l-2a)(g_L- s ± ) 2 



(2tt) 2 (2tt) 2 Jo ql [p + a (i _ a )q 2 x + L 2 {m 2 )} 



[(37) 
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FIG. 4: Feynman diagram at the first non-trivial order that could potentially contribute to the twist-3 quark-gluon correlation 
functions. 



which also vanishes from the symmetry of the da integration. We thus verify that p(x, x) = 0. 

To conclude this subsection on the calculation with a scalar diquark, we make a comment on the contribution to 
the quark-gluon correlation function from the diagram in Fig. 2J In order to get the SSAs, as discussed earlier, we 
need a spin flip between the two partonic states on the opposite side of the cut in the diagram in Fig. 0] Since the 
quark-gluon composite state on the left was initiated from a single quark state, the spin flip contribution can only 
come from the mass term of the quark. Therefore, the contribution of the diagram in Fig. [4] to the quark-gluon 
correlation functions relevant to the SSAs is expected to be proportional to the quark mass and is therefore small. 
Our explicit calculation shows that the diagonal correlation function T q ^p(x,x) relevant to the leading gluonic pole 
contribution vanishes at this order which is consistent with the fact that when the gluon momentum vanishes, there 
is no spin flip between two quarks on the opposite side of the cut. For the off-diagonal correlation functions relevant 
to the fermionic pole contribution, we find 



T (s j(0,x) . =-T« (0,aO . m = NcC ™ K {l-x) m {m + xM) 2 



-.(») 



Fig.H 



Fig. g] 87T 

d 2 k ± d 2 q ± F(k 2 ) 2 



(2tt) 2 (2tt) 2 [ q l + m 2] [ fc 2 + L 2 {m 2 f 



(38) 



which is clearly proportional to the mass of quark. 



C. Calculation with an axial-vector diquark 



In order to test the sensitivity of our results derived in the last subsection on the choice of the scalar diquark, we 
present in this subsection quark-gluon correlation functions calculated with an axial-vector diquark. 

By using the same Feynman diagram in Fig. [3] with the Feynman rule for an axial- vector spectator, and the same 
cut vertices, we derive the quark-gluon correlation functions relevant to both leading gluonic and fermionic pole 
contribution to the SSAs. Since the pole structure of the Feynman diagram in Fig. [3] is insensitive to whether the 
spectator is a scalar or an axial-vector, we find, like the case of a scalar diquark, that all off-diagonal quark-gluon 
correlation functions relevant to the leading fermionic pole contribution vanish, 

T$(0,x)=T£l F (0,x)=0. (39) 

For the diagonal quark-gluon correlation functions relevant to the leading gluonic pole contribution, we obtain 

T W (x x) N c C F gX 2 v g v f d*k ± rfV [gj - ( g± ■ s^) 2 } F (k 2 )F((k + g) 2 ) 

' j " in X[m + XM) J (27T) 2 (2n) 2 q l[kl+L 2 (m 2 )}[(k±+q±) 2 +L 2 (m 2 )}' [W) 

which is the same as that in Eq. (f2"U)) except the overall (1 — x) factor is replaced by x due to the difference in diquark 
spin. Therefore, the rest of derivation and discussion in the last subsection following Eq. (|2"9"|) should be the same for 
the case of an axial-vector diquark. We find 



T >^ n - ^^<rn + X M) [ |% 1 ln ^±j^) + ^ , (41) 

«' point-like 16tH J (2tt) 2 q\ + AL 2 s (m 2 ) y/ q\ + AL 2 (m 2 ) - qj_ 

0">U, = ^1^(1-*)V + *M) (jjIj) 2 , (42, 
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which are the same as those in Eqs. (l3Tj) and (|33|) except that one factor of (1 — x) is replaced by x. We also explicitly 
verify that T AqF (x,x) = when it is calculated with an axial-vector diquark. 

To conclude this section, we summarize our key results as follows. We find, in terms of an explicit calculation in 
the quark-diquark model of the nucleon, that at the leading non-trivial order all quark-gluon correlation functions 
relevant to the leading fermionic pole contribution to the SSAs vanish, 

T qiF (0,x) = T q>F {x,0) = 0, T Aq M0,x)^-T AqtF (x,0)=0. (43) 

We also verify that T AqF (x,x) = 0, and find that only the diagonal quark-gluon correlation function, T qF (x,x), is 
finite. 



IV. CONNECTION TO TMD PARTON DISTRIBUTION FUNCTIONS 



As we stressed in the introduction of this paper, the collinear factorization approach and the TMD factorization 
approach to the SSAs are closely connected and complementary to each other. It was shown in terms of their 
operator definitions that the twist-3 quark-gluon correlation function T q ^ F (x,x) is related to the moment of quark 
Sivers function fxr(x,kj_) 0, 

T q A*>*) = ]jjf / d 2 k ± klft T (x,kl), (44) 

where f^p(x, k\) is the quark Sivers function defined via the Drell-Yan process, which is related to the quark Sivers 
function defined in the semi- inclusive deep inelastic scattering by a minus sign [17l . [l8j . In this section, we explicitly 
verify this relation in Eq. (|44[) by comparing the quark-gluon correlation function T q>F (x,x) calculated in this paper 
with the Sivers function calculated in the same quark-diquark model of the nucleon. 



q 




FIG. 5: Lowest order Feynman diagram for the quark Sivers function in the quark-diquark model of the nucleon [41 



The quark Sivers function in the quark-diquark model of the nucleon has been calculated in Ref. [41] . The Feynman 
diagram to the lowest non-trivial order for the quark Sivers function is shown in Fig. O In terms of an explicit 
calculation with a scalar diquark we obtain 



At \ x > ^x) 



N c C F g\ 2 s g s (1 - x)M(m + xM) ^ k\ + L 2 {m 2 ) 



poi„t-iikc 4(2tt) 4 kl[kl+L 2 s (m 2 )} "* L 2 (m 2 

for a point-like interaction between the nucleon, the quark, and the spectator diquark, and 



/it ( x i 



dipolar 



N c C F g\ 2 s g s 
4(2tt) 4 



(1 - x) 3 M(m + xM) 



[A 



212 



L 2 (A 2 )[kl + L 2 (A 2 )] 



(45) 



(46) 



for using the dipolar form factor for the interaction between the nucleon, the quark, and the spectator diquark. 
Similarly, we find, in terms of the calculation with an axial- vector diquark, 



/it 



N c C F g\ 2 v g v xM(m + xM) ^k 2 ± + L 2 s (m 2 ) 



point — like 



4(2vr) 4 k\ [k\ +L 2 s (m 2 )} 



and 



/it ( x ik±) 



dipolar 



N c C F g\lg v 2 

— — x(l - x) M (m + xM) 



L 2 (m 2 ) 

[A-s] 2 



L 2 (A 2 )[kl+L 2 (A 2 )T 



(47) 



(48) 
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respectively. The calculated results here are the same as those obtained in Ref. [41| except the overall constant factor 
[A 2 ] 2 for those with the dipolar form factor. The difference is, as explained earlier, due to a slightly different choice 
of the form factor so that the calculated twist-3 correlation functions as well as the quark Sivers function appear to 
have the right dimension. 

In order to verify the relation in Eq. (|44|) . we need to take the moment of the quark Sivers functions calculated above 
in the same quark-diquark model of the nucleon. However, the moment of the quark Sivers functions in Eqs. (|45[) and 
(1471) calculated by using the point-like interaction is logarithmically divergent, for example, 



1 

M 



point— like 



N c C F g\ 2 s g s 
16tt 2 
d 2 kt 



(1 - x)(m + xM) 



(2tt) 2 [k\ + L 2 s {m 2 )\ 



In 



k\ + L 2 (m 2 ) 
Ll{m?) 



(49) 



The moment of the quark Sivers function in Eq. (|49j) is clearly not necessary to be the same as the twist-3 quark- 
gluon correlation function in Eq. (1311) , even if one imposes the same ultraviolet cutoff on the transverse momentum 
integration in both Eqs. (|3"Tj) and (|49[) . This is because the ultraviolet divergence of the calculated quark-gluon 
correlation function in Eq. (|31[) was not regularized in the same way as that in the Sivers function calculation in 
Eq. (1451) . This example indicates that if one wants to compare the both sides of Eq. (1441) perturbatively by projecting 
the equation onto a parton state, one has to specify the regularization and renormalization condition for the ultraviolet 
divergence in both sides. In general, the relation in Eq. (|44[) is not necessarily valid for all orders in perturbative 
calculations if one does not regularize and renormalize the ultraviolet divergence in the same way for the both sides. 

If we regularize and renormalize the ultraviolet divergence in both sides of Eq. (144]) in the same way, we should 
expect the relation to hold. To explicitly demonstrate this, we compare the quark-gluon correlation functions and the 
quark Sivers function calculated with the same dipolar form factor, 



scalar diquark: 



axial- vector diquark: 



1 

17 
i 

17 



d 2 k±k'j_f 1 rp (x, fcj_) 



d 2 k±k 2 _f 1 r^ \x, k\) 



N c C F g\ 2 s g s 
16(2tt) 3 

N c C F g\ 2 v g v 
16(2tt) 3 



(1 - xf(m + xM) 



A 2 



x(l - x) 2 (m + xM) 



A 2 



(50) 



(51) 



The right-hand-side of above equations are clearly equal to the quark-gluon correlation functions in Eqs. (|33[) 
(H2J), respectively. 



and 



V. SUMMARY AND CONCLUSIONS 



In this paper, we calculate various twist-3 quark-gluon correlation functions of a transversely polarized nucleon in 
the quark-diquark model of the nucleon. Our calculations are done with the diquark being a scalar particle as well 
as being an axial-vector particle. We have found from our calculation at the first non-trivial order that all quark- 
gluon correlation functions relevant to the leading fermionic pole contribution, T qj p(0,x), T^ qt p(0,x), T q ^p(x,0), and 
Ia 9 ,f(^,0), vanish. Only one of the diagonal quark-gluon correlation functions relevant to the leading gluonic pole 
contribution, T qi p(x, x), is finite. The other diagonal quark-gluon correlation function, T& qt p(x, x), also vanishes from 
both the symmetry argument and explicit calculation. Our conclusions are independent of the diquark being a scalar 
or an axial- vector. 

Although our results are derived from a specific model calculation, the features of the calculated results should 
allow us to conclude with confidence that the diagonal quark-gluon correlation function T q j?(x, x) is much larger than 
all other quark-gluon correlation functions that are relevant to the leading soft pole contribution to the SSAs. This 
conclusion is significant and important for phenomenological study of the SSAs. It enables us to study the physics 
of SSAs without including too many unknown correlation functions at the early stage of probing this new domain 
of QCD dynamics. However, it is the limitation of the quark-diquark model of the nucleon that we are not able to 
calculate the tri-gluon correlation functions in this model, which are closely connected to the quark-gluon correlation 
functions via perturbative evolution [33l l36j . 

As we explained in last section, it requires a caution in using the relation between the twist-3 quark-gluon correlation 
function T qt p(x,x) and the moment of the quark Sivers function in Eq. (|44p. Since both sides of the equation, the 
twist-3 quark-gluon correlation function T q ^p{x,x) on the left and the moment of quark Sivers function on the right, 
are perturbatively divergent, the relation is valid only if the same regularization and renormalization scheme is adopted 
to the calculation of both sides. As an example, we demonstrate in the last section that the relation could be violated 
perturbatively if different regularization and renormalization schemes were used; and the relation is valid if the same 
scheme were used in both sides. 
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